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ABSTRACT 

Closely  spaced  observations  of  nonlinear  internal  waves  (NLIWs)  were  made  on  the  outer  continental  shelf 
off  New  Jersey  in  June  2(X)9.  Nearly  full  water  column  measurements  of  current  velocity  were  made  with  four 
acoustic  Doppler  current  profilers  (ADCPs)  that  were  moored  about  5  km  apart  on  the  bottom  along  a  line 
approximately  normal  to  the  baihymetry  between  water  depths  of  67  and  92  m.  Density  profiles  were  obtained 
from  two  vertical  strings  of  temperature  and  conductivity  sensors  that  were  deployed  near  each  of  the  interior 
ADCP  moorings.  In  addition,  a  towed  ScanF  ish  provided  profiles  and  fixed-level  records  of  temperature  and 
salinity  through  several  NL1 W  packets  near  the  moorings.  Several  case  studies  were  selected  to  describe  the 
propagation  of  lhe  NLIWs.  One  to  three  solitary  waves  of  depression  were  observed  in  five  selected  packets. 
There  were  also  occurrences  of  muliiple-phase  dispersive  wave  packets.  The  average  propagation  speed 
corrected  for  adveetion  of  the  observed  waves  was  0.51  ±  0.09  ms1.  The  waves  were  directed  primarily 
shoreward  (~- northwestward)  along  the  mooring  line  with  average  wavelengths  and  periods  of  about 
300  ni  and  10  min.  respectively.  Wave  amplitudes  and  energies  decreased  with  decreasing  water  depth.  The 
observed  wave  parameters  tan  be  locally  described  by  a  two-layer  Korteweg-dc  Vries  <KdV)  model,  except 
for  the  decreasing  amplitudes,  which  may  be  due  to  shear-induced  dissipation  and/or  bottom  drag.  The 
various  complementary  observations  utilized  in  this  study  present  a  unique  description  of  NLIWs. 


1.  Introduction 

Nonlinear  internal  wave  (NL1W)  packets  arc  com¬ 
monly  observed  on  continental  shelves  throughout  the 
world.  They  are  commonly  thought  to  evolve  from  in¬ 
ternal  tides  generated  by  the  interaction  of  the  baro- 
tropic  tide  with  shelf-edge  topography  (Tang  ct  al.  2007; 
Colosi  et  al.  2001)  and  to  propagate  shoreward  when  the 
interaction  is  with  the  shelfbreak.  Characteristics  of  the 
waves,  which  are  often  described  as  packets  of  solitons, 
change  as  the  water  depth  decreases.  Their  amplitudes 
can  be  quite  large,  and  the  waves  can  significantly  im¬ 
pact  acoustic,  optical,  and  biogeochemical  properties  of 
the  water  column.  A  variety  of  techniques  have  been 
used  to  observe  NLIWs  such  as  using  moored  temper¬ 
ature  and  conductivity  sensors,  current  meters  (Colosi 
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et  al.  2001;  Duda  et  al.  2004;  Hallock  and  Field  2005; 
MacKinnon  and  Gregg  2003),  and  high-frequency  pres¬ 
sure  sensors  (Mourn  and  Nash  2008).  Shipboard  mea¬ 
surements  utilizing  acoustic  Doppler  current  profilers 
(ADCPs),  conductivity-temperature-depth  (CTD)  pro¬ 
filers,  and  echo  sounders  have  been  used  to  make  multi 
pie  transects  across  wave  packets  (e.g.,  Shroyer  et  al.  2008). 
NLIWs  have  been  particularly  well  documented  through 
satellite  synthetic  aperture  radar  (SAR)  imagery  (Jackson 
and  Apcl  2(X)4).  NLIWs  can  have  a  surface  signature  de¬ 
tectable  by  both  ship  and  satellite  radars  in  low  wind 
conditions. 

NLIWs  were  intensely  studied  as  part  of  the  Shallow 
Water  Acoustics  Experiment  in  the  summer  of  2006 
(SW06)  (Tang  et  al.  2007).  This  large  multidisciplinary, 
multi-institution  project  off  the  coast  of  New  Jersey  fo¬ 
cused  on  a  hierarchy  of  processes,  including  internal  tides, 
jets,  density  fronts,  and  NLIWs,  that  could  affect  U.S. 
Navy  acoustic  systems.  SW06  objectives  are  focused  on 
the  structure,  intensity,  and  direction  of  the  internal  wave 
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FIG.  1.  (a)  The  location  of  the  experi- 
menl  is  marked  by  the  box.  (b)  Detail 
showing  mooring  locations:  ADCPs  A I 
A4  are  marked  by  triangles:  string  moor¬ 
ings  SI  and  S2  are  collocated  with  ADCPs 
A2  and  A3,  respectively:  and  relevant 
ScanFish  lows  are  indicated  by  track  lines 
about  2  km  northeast  of  the  mooring  line, 
(c)  Bottom  profile  along  the  mooring 
line:  the  ADCP  locations  are  indicated. 


Distance  (km) 


energy  flux  at  the  shelf  break.  During  SW06,  62  moor¬ 
ings  were  deployed  in  a  “T"  geometry,  and  5  of  these 
moorings  contained  ADCPs,  which  provided  nearly  full 
water  column  velocity  profiles.  Three  of  these  ADCP 
moorings  were  placed  about  20  km  apart  along  the  likely 
propagation  direction  of  the  NLlWs  while  three  were 
10—15  km  apart  along  a  line  perpendicular  to  the  prop¬ 
agation  direction.  In  addition,  multiple  ships  were  used 
to  chase  the  packets,  and  gliders  provided  background 
observations.  Shipboard  observations  included  ADCP, 
X-band  radar,  temperature,  salinity,  and  turbulence  mea¬ 
surements. 

Analyses  of  satellite  imagery'  suggest  that  the  width 
of  a  typical  NLIW  packet  can  range  from  several  to 
more  than  5  km  (Porter  and  Thompson  1999).  Thus  far. 
this  resolution  has  only  been  approached  using  a  ship 
in  chase  mode  collecting  ADCP  and  CTD  data.  Moored 
instruments  spaced  at  about  5  km  or  less  would  allow 
coherent  measurements  of  waveforms  and  direct  deter¬ 
mination  of  phase  and  group  velocities.  Individual  wave 
packets  could  be  analyzed  in  detail  for  modal  content 
and  propagation  characteristics.  As  part  of  the  Naval 
Research  Laboratory's  (NRL's)  field  exercise  near  the 
shelf  break  off  New  Jersey,  four  ADCP  moorings  were 
spaced  about  5  km  apart  along  a  line  approximately 


normal  to  the  bathymetry  (Fig.  1).  The  field  measure¬ 
ment  were  performed  as  part  of  NRL  Magnetic  Ocean 
Forecasts  project,  whose  main  objective  was  to  generate 
estimates  of  the  ocean  contribution  to  the  magnetic  noise 
from  sources  such  as  NLIWs.  This  topic  is  not  addressed 
here.  However,  the  field  measurements  provide  the  op¬ 
portunity  to  analyze  closely  spaced  mooring  data  and  to 
fill  in  gaps  of  our  understanding  of  NLIW  dynamics  not 
possible  from  earlier  field  studies.  This  paper  focuses  on 
analyses  of  several  packets  of  large  NLIWs  that  propa¬ 
gated  through  our  mooring  array. 

In  section  2,  w'e  describe  the  experiment  and  the  ob¬ 
servations  made.  1  n  section  3,  analyses  of  selected  events 
in  the  data  records  and  determination  of  wave  charac¬ 
teristics  are  presented.  Section  4  results  arc  interpreted 
using  first-order  Kortewcg-de  Vries  (KdV)  theory  for 
a  two-layer  system.  Some  conclusions  are  presented  in 
section  5. 

2.  Observations 

An  extensive  set  of  oceanographic  measurements  was 
made  by  the  Naval  Research  Laboratory  between  18  June 
and  1  July  2009  (corresponding  to  yeardays  168  and  181) 
during  an  R/V  Occanus  cruise.  The  study  area  was 
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centered  at  about  39.5°N  and  72.7°W,  about  75  miles  cast 
of  New  Jersey  (Fig.  la).  The  instrumentation  consisted 
of  four  ADCP  m(X>rings;  two  string  moorings  with  tem¬ 
perature  and  conductivity;  and  a  towed  undulating  body, 
SeanFish.  Eight  magnetometers  were  also  deployed  on 
the  bottom  next  to  the  ADCP  moorings  but  are  not  dis¬ 
cussed  in  this  paper.  After  all  of  the  moorings  were  de¬ 
ployed,  ship  observations  were  made  along  a  line  parallel 
to  the  mooring  line,  offset  about  2  km  toward  the 
northeast  (Fig.  lb).  The  mooring  line  was  avoided  in  or¬ 
der  to  not  disturb  the  magnetic  measurements. 

a .  ADCP  moorings 

Four  ADCP  moorings  were  deployed  on  the  outer 
continental  shelf  off  New  Jersey  and  southwest  of  the 
Hudson  Canyon  at  water  depths  ranging  between  64  and 
92  m  (Fig.  1 ).  The  bottom  profile  along  the  mooring  line 
is  shown  in  Fig.  le.  Mooring  A1  is  located  at  a  fairly  flat 
portion  of  the  bottom,  A2  is  located  at  the  edge  of  the 
Hat  portion  where  the  bottom  profile  begins  to  rapidly 
deepen,  and  A3  and  A4  are  located  on  the  sloping  por¬ 
tion  of  the  profile.  The  moorings  were  deployed  for 
about  9  days  from  19  June  2(X)9  (yearday  169)  along  a 
line  oriented  at  310°  relative  to  true  north  and  were 
approximately  5  km  apart.  The  mooring  line  is  approxi¬ 
mately  normal  to  the  bathymetry.  Numerous  SAR  images 
collected  in  this  region  during  August  and  September 
2006  in  conjunction  with  SW06  suggest  that  NLIW 
crests  arc  likely  parallel  to  the  isobaths  along  our  moor¬ 
ing  line. 

The  ADCPs  were  307.2-kHz  broadband  “Workhorse 
Sentinels"  manufactured  by  RD  Instruments.  They 
were  deployed  in  trawl  resistant  mounts  called  “BarnysM 
after  their  barnaele-like  shape  (Perkins  et  al.  2(XK1)  that 
were  lowered  rather  than  free  dropped  to  the  bottom. 
They  were  then  positioned  so  that  the  tilt  of  the  trans¬ 
ducer  head  from  the  vertical  (ideal)  orientation  was 
minimized.  The  moorings  were  placed  at  depths  ranging 
from  about  67  to  92  m,  with  the  seaward  end  of  the 
mooring  line  near  the  shelf  break  of  about  100  m.  The 
ADCP  heads  were  situated  about  0.5  m  off  the  bottom 
and  recorded  current  profiles  from  about  4  m  off  the 
bottom  to  within  about  4  m  of  the  surface.  Velocity 
measurements  were  recorded  at  a  vertical  resolution  of 
1  m  every  minute  using  a  sampling  scheme  that  used  70 
pings  distributed  over  the  sampling  period.  The  accu¬ 
racy  is  0.5%  ±  0.5  cm  s  1  and  the  random  error,  w  hieh  is 
dependent  on  the  depth-bin  size,  is  reported  as  a  stan¬ 
dard  deviation  of  1.6  cm  s-1.  Velocities  were  corrected 
for  local  magnetic  variation.  Each  mooring  also  con¬ 
tained  sensors  for  measuring  temperature  and  pressure. 
The  ADCP  velocities  are  rotated  (39.5°  clockwise)  for 
the  analyses  here  so  that  the  east-west  (U)  velocities  are 


oriented  along  the  mooring  line,  where  negative  U  is 
onshore  and  positive  U  is  offshore. 

h.  Siring  moorings 

Two  temperature-conductivity  sensor  strings  (SI  and 
S2)  were  anchored  on  the  bottom  beneath  a  surface 
buoy  approximately  300  m  from  the  two  interior  ADCP 
moorings  (A2  and  A3;  Fig.  1).  Each  string  contained 
10  Sea-Bird  Electronics  (SBE)  37  serial  interface 
MieroCATs,  high-accuracy  conductivity,  and  tempera¬ 
ture  sensors.  About  half  of  the  MieroCATs  on  each 
string  contained  pressure  sensors.  They  were  distributed 
in  a  telescoping  scheme,  starting  with  a  spacing  of  3  m  at 
about  8  m  below  the  surface  and  ranging  to  a  spacing  of 
10  m  al  about  10  m  above  the  bottom.  Temperature  is 
accurate  to  about  0.002°C,  and  conductivity  is  accurate 
to  about  0.003  mS  cm 

c.  SeanFish 

Repeated  tows  of  a  SeanFish  MK  II  manufactured  by 
EIVA  (http://www.eiva.dk)  were  made  along  a  line  par¬ 
allel  to  the  mooring  line  (Fig.  1)  throughout  the  cruise 
period,  shifted  by  2  km  northeast  of  the  mooring  line  to 
minimize  interference  with  the  deployed  magnetome¬ 
ters.  The  SeanFish  is  a  towed  undulating  vehicle  that 
can  fly  in  undulating,  fixed  depth,  and  terrain  following 
modes.  It  records  numerous  physical  and  optical  pa¬ 
rameters,  but  just  the  temperature,  conductivity,  and 
pressure  data  are  used  in  analyses  here.  The  SeanFish 
was  typically  towed  at  5.5  kt  throughout  the  cruise.  Both 
level  and  undulating  modes  were  used  to  sample  NLIW 
packets. 

d.  Hydrography 

CTD  stations  were  taken  at  mooring  sites  A1-A3  on 
20  June  (yearday  170)  and  near  mooring  sites  A1-A4 
along  the  SeanFish  line,  approximately  2  km  northeast 
of  the  moorings,  on  25  June  (yearday  175).  Each  CTD 
station  consisted  of  three  consecutive  downcasts,  made 
to  just  above  the  bottom.  On  20  June,  the  mixed  layer 
was  shallow,  5-10  m  in  depth,  and  the  seasonal  ther- 
mocline  ranged  from  about  10  to  25  m  in  depth.  Well- 
mixed  bottom  layers  that  were  about  5  m  thick  were 
observed  at  Al  and  A2.  At  A3,  a  well-mixed  layer  aboul 
15  m  thick  was  located  just  above  a  near-bottom  density 
gradient  where  density  <tx  changed  by  0.3  kg  m  *  over 
the  bottom  5  m.  Brunt-Vaisalii  frequencies  ranged  from 
26  to  31  eph  in  the  seasonal  thermoeline  and  from  16  to 
19  eph  near  the  bottom.  By  25  June,  the  mixed  layer  had 
eroded  and  the  seasonal  thermoeline  weakened  and 
deepened.  Well-mixed  bottom  layers  disappeared  at  A4, 
A3,  and  A2.  However,  the  bottom  layer  at  Al  thickened 
to  about  15  m  and  was  less  dense  by  about  0.2  kg  m 
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A1 


A2 


A3 


A4‘ 


PlG.  2.  Average  profiles  of  <r(  based  on  CTD  easts  conducted  at  (or  near)  ADCP  sites  al  deployment  lime,  t  he 
panel  A4'  is  so  labeled  because  il  was  located  about  2  km  northeast  of  A4.  Horizontal  dotted  lines  indicate  the  depth 
of  maximum  variance  of  vertical  velocity  from  ADCP  records  and  hence  the  mean  depth  of  the  upper  layer  in  a 
2-layer  system.  Vertical  lines  indicale  vertical  layer  averages  of  density.  The  red  and  green  curves  are  2.4-h  case  111 
and  IV  average  profiles,  respectively,  of  from  string  moorings;  diamonds  indicate  instrument  locations. 


from  the  profiles  taken  5  days  earlier.  Similarly,  the 
density  decreased  near  the  bottom  at  the  other  moor¬ 
ings.  Brunt-VaisaJa  frequencies  were  much  lower,  and 
there  were  no  well-defined  peaks  for  the  profiles  taken 
on  25  June.  In  all  cases,  variability  was  highest  over  the 
three  consecutive  casts  in  the  seasonal  thermoelinc  and 
is  indicative  of  internal  wave  activity.  Results  of  the 
analyses  below  are  compared  to  a  two  layer  dynamical 
system.  We  determine  suitable  two-laver  parameters 
using  density  profiles  from  the  described  CTD  easts. 
Upper-layer  depth  hx  is  taken  as  the  depth  of  maximum 
time  variance  of  vertical  velocity  (W)  at  each  ADCP 
site;  h2  is  then  the  difference  between  the  water  depth 
//  and  h\.  At  Al,  A2,  A3,  and  (near)  A4,  three  succes¬ 
sive  CTD  casts  were  conducted  yielding  an  average 
density  profile  (Fig.  2)  at  each  site.  Layer  densities  are 
taken  as  the  respective  vertical  averages  of  density 
within  the  layers;  these  averages  arc  indicated  in  Fig.  2. 
Superimposed  on  the  profiles  at  A2  and  A3  are  mean 
profiles  from  the  moored  CTD  strings  (SI,  S2)  for  2.4-h 
periods  on  21  and  26  June  (cases  III  and  IV;  see  Table 
1).  The  latter  of  these  (case  IV)  further  shows  the 
temporal  change  in  stratification  over  the  course  of  the 
experiment. 


3,  Analysis 

During  the  2-week  cruise,  six  cases  of  strong  NLlWs 
were  identified.  Three  were  identified  at  the  beginning 
of  the  cruise,  and  three  were  identified  near  the  end  of 
the  cruise.  The  middle  of  the  cruise  period  was  at  the 
peak  of  the  spring  tide  and  no  strong  NLIW  activity  was 
observed;  the  length  of  the  experiment,  however,  is  too 
short  to  make  an  inference  on  this  issue.  The  six  case 
studies  arc  summarized  in  fable  1  according  to  the 
type  of  measurement  made.  The  string  moorings  were 
not  in  place  for  cases  I  and  II,  and  all  moorings  had 

Table  I.  Periods  of  energetic  NLIW  activity  were  measured  by 
ADCP  moorings  A I  A4.  temperature -conductivity  string  moor- 


ings  SI  and  S2,  and  level  ScanEish  tows.  The  day  indicates  decimal 
day.  The  number  of  events  observed  at  each  location  is  shown. 

C  asc/day 

Al 

A2 

A3 

A4 

SI  S2 

Level  tow 

1/ 1 70.3 

2 

2 

2 

11/170.7 

3 

3 

3 

3 

tit/171.3 

2 

2 

2 

2 

2 

IV/I77.0 

2 

2 

2 

2 

2 

V/ 177.9 
VI/180.5 

I 

1 

1 

l 

Multiple 
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Flo.  3.  Vertical  velocily  H'ai  23  m  ai  ADCP  2.  Vertical  dolled  lines 
indicale  lhe  5  cases  for  lhe  ADCP  soliton  evenls. 


been  recovered  by  the  time  cases  V  and  VI  were  ob¬ 
served.  Packets  were  not  evident  at  A4  during  cases  I 
and  IV, 

Vertical  velocity  from  ADCP  2  at  the  depth  of  max¬ 
imum  variance  is  plotted  in  Fig.  3,  where  a  number  of 
large-amplitude  events  are  evident  in  the  record.  A  few 
of  these  were  selected  for  analysis  based  on  three  cri¬ 
teria:  relatively  large  amplitude;  resemblance  to  simple 
internal  solitary  waves;  and  appearance  in  at  least  three 
of  the  ADCP  records,  thus  allowing  tracking  and  prop¬ 
agation  speed  determinations.  The  mooring  line  was 
chosen  to  be  parallel  to  expected  NLIW  propagation 
direction.  Clearly,  this  direction  can  vary  for  different 
waves  or  packets  of  waves  passing  a  given  point.  If  we 
take  the  .v  along  the  mooring  line  and  0  as  the  departure 
of  wave  direction  from  the  x  axis,  then  (for  plane  wave 
fronts)  the  apparent  phase  speed  (along  .v)  C'  =  C/cos0, 
where  C  is  the  actual  phase  speed.  This  means  that  the 
observed  speeds  will  be  greater  than  or  equal  to  the 
actual  speed.  For  0  up  to  about  25°,  the  ratio  CfC  is 
greater  than  0.9,  so  small  directional  uncertainties  will 
not  result  in  large  errors  in  estimates  of  C.  Five  cases 
(1-V;  see  Table  I),  each  consisting  of  1-3  closely  spaced 
solitary  wave  events,  are  indicated  in  Fig.  3.  A  total  of 
10  solilons  were  tracked  through  most  of  the  moorings, 
resulting  in  36  observations  at  the  four  ADCP  locations 
and  4  observations  at  SI  and  S2.  Case  I  and  IV  solitons 
were  not  seen  at  A4,  The  case  IV  solitons  (2)  were  also 
observed  twice  with  the  ScanFish.  The  prominent  event 
near  yearday  172.5  is  predominantly  a  packet  of  large- 
amplitude  internal  waves  rather  than  well-defined  soli¬ 
tary  waves.  Case  VI  consisted  of  two  ScanFish  tows,  near 
the  A2  location,  conducted  after  recovery  of  all  moor¬ 
ings.  We  analyze  the  observations  in  the  light  of  a  two- 
laver,  first-order  KdV  model  (Small  2001;  Holloway 
1987;  Gear  and  Grimshaw  1983). 


FlO.  4,  Schemalicof  KdV  solilon  of  depression  (red)  and  associated 
vertical  velocily  (black). 


a.  Phase  speeds 

A  number  of  factors  determine  apparent  propagation 
speeds,  along  a  particular  direction,  of  NLIWs  (Small 
2001;  Gear  and  Grimshaw  1983;  Phillips  1969).  These 
include  stratification,  water  depth,  wave  amplitude,  ac¬ 
tual  propagation  direction,  and  advection  (Doppler  shift) 
and  shear  by  the  ambient  flow  in  which  the  waves  are 
embedded. 

Phase  speeds  of  a  number  of  selected  NLIWs  ob¬ 
served  in  the  ADCP  velocities  are  estimated  using  the 
times  that  a  particular  event  passes  each  ADCP  and  the 
distance  between  the  moorings.  Because  the  internal 
waves  are  superimposed  on  the  ambient,  depth-averaged 
low-frequency  flow',  mostly  the  barotropic  tide,  a  Doppler 
correction  is  made  so  that  calculated  speeds  are  relative 
to  this  ambient  flow.  The  events  analyzed  were  selected 
by  examining  records  of  vertical  velocity  W  that  were 
robust  and  appeared  in  at  least  two  data  records.  The 
characteristic  shape  of  a  first-order  KdV  solilon  of  de¬ 
pression  (Holloway  1987;  Gear  and  Grimshaw  1983)  is 
described  by 

£{xj)~a  scch2(^—£—\  and 


d£  2a  C  2  (x 

W  -  —  ~  — —  sech 
dt  L 


(1) 


where  £  is  elevation  (interface  displacement  for  a  two- 
layer  system),  a  <  0  for  a  soliton  of  depression,  L  is 
a  length  scale,  and  C  is  propagation  speed.  For  first- 
vertical -mode  waves,  the  horizontal  velocily  is  propor¬ 
tional  to  £  and  changes  sign  near  the  depth  of  maximum 
W.  Hence,  U can  also  be  described  bv  the  sech2  function. 
For  Eulerian  measurements,  we  can  set  x  -  0  and  define 
a  time  scale  r  =  L/C  and  reference  time  /0  to  simplify  the 
following.  These  functions  are  plotted  in  Fig.  4  with  /,>  = 
50  min  and  r  =  10  min.  The  period  of  the  solilon  can  be 
defined  as  3.6r,  the  point  at  which  the  magnitude  of  £ 
falls  to  10%  of  the  maximum  (Holloway  1987);  likewise. 
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FlG.  5.  Vertical  vek>cily  from  ADCP  A1  (case  11).  The  abscissa  is 
in  minutes  for  easy  comparison  with  the  analytical  plot  in  Fig.  4 
Diamonds  on  lhe  curve  denote  actual  data  samples  (1-min  interval). 


a  wavelength  is  defined  as  A  =  3.6L.  It  should  he  noted 
that  the  amplitude  of  W  is  proportional  to  ah.  When 
a  soli  ton  of  depression  arrives  at  a  moored  measuring 
system,  W  first  goes  negative  to  a  minimum  value,  then 
crosses  zero  to  a  maximum,  and  subsequently  decays 
back  to  zero.  We  define  the  time  of  arrival  of  a  soliton 
event  at  a  mooring  as  (0,  where  IV  crosses  zero  from  the 
minimum  to  the  maximum.  An  example  of  such  an  event 
in  a  moored  current  record  is  shown  in  Fig.  5.  Time  in 
relative  minutes  is  shown  to  more  easily  compare  with 
Fig.  4.  The  actual  value  of  (q  is  yearday  170.85. 

T  he  event  in  Fig.  5  was  detected  earlier  in  the  other 
three  ADCP  records,  from  locations  southeastward 
from  Al.  There  were  also  two  additional  events  fol¬ 
lowing  the  one  shown  that  appeared  in  all  four  records. 
The  northwestward  progression  of  these  events  is  shown 
in  Fig.  6;  a  somewhat  longer  time  interval  for  all  four 
records  are  plotted  to  include  the  three  events  (labeled 
a-c).  The  times  of  the  events  and  distances  between  the 
moorings  (—5  km)  yield  apparent  propagation  speeds 
between  respective  mooring  locations.  The  three  events 
moving  from  A4  to  A3,  from  A3  to  A2.  and  from  A2  to 
Al  yield  nine  estimates  of  propagation  speed.  These 
range  from  -0.74  to  -0.44  m  s  '.The  calculated  speeds 
are  interpreted  as  soliton  phase  speed  C  superimposed 
on  the  ambient  flow  (due  primarily  to  the  barotropic 
tide).  The  ambient  flow'  is  estimated  by  low-pass-filtered, 
depth-averaged  currents  (LF  currents)  from  the  ADCP 
records.  The  half-amplitude  point  of  the  filter  is  set  at  a 
period  of  2  h.  For  illustration,  a  further  average  of  ro¬ 
tated  LF  current,  over  the  four  ADCP  records,  is  plotted 
in  Fig.  7.  Fluctuations  are  due  mainly  to  the  M2  tide,  but 
there  is  also  a  southeastward  mean  as  well.  Actual  cor¬ 
rections  for  speed  calculations  were  made  as  follows,  for 
example:  An  event  passes  from  A3  to  A2  in  a  time  in¬ 
terval  A/  over  the  distance  Aa*  (~5  km).  We  calculate  an 
average  LF  current  t/32,  where  the  double  overbar  in¬ 
dicates  averaging  the  LF  flow'  over  A2  and  A3  and  over 


A/.  Then,  we  have  Aar  =  (C  +  )At  Solving  for  C yields 
a  corrected  phase  speed.  For  the  events  in  Fig.  6,  C 
ranges  from  -0.65  to  -0.48  ms  1 .  Here,  C  is  an  estimate 
of  the  intrinsic  internal  wave  phase  speed.  We  note  that 
this  correction  docs  not  address  mean  vertical  shear, 
which  can  also  affects  the  phase  speeds  and  other  wave 
properties  (Alford  et  al.  2010). 

In  all,  five  cases  (of  which  the  above  is  case  II)  cap¬ 
tured  by  the  ADCP  moorings,  some  with  subcases  (a,  b. 
etc.),  were  selected  for  a  total  of  10  solitons  passing 
at  least  three  of  the  moorings.  Each  of  the  36  events 
observed  resemble  Fig.  5  and  are  hence  solitons  of  de¬ 
pression.  Means  and  standard  deviations  (for  each  lo¬ 
cation)  of  both  advcction-correctcd  and  uncorrected  C 
estimates  are  plotted  in  Fig.  8.  The  uncorrected  speeds 
show  an  increase  in  (negative)  amplitude  with  increas¬ 
ing  depth;  the  variability  suggests  that  this  may  not  be 
a  significant  change.  The  corrected  speeds,  however, 
show  less  change  and  do  not  follow  the  depth  profile.  We 
conclude  that  the  changes  arc  not  significant  and  that 
the  composite  mean  of  -0.51  ±  0.09  m  s  1  is  a  best  es¬ 
timate  of  the  phase  speed  along  the  mooring  line,  based 
on  the  selected  events  from  the  ADCP  observalions. 
Furthermore,  this  mean  value  is  not  significantly  differ¬ 
ent  from  the  mean  for  the  uneorrected  speeds  (-0.54  ± 
0.11  m  s-1).  Horizontal  current  velocities  associated 
with  the  wave  events  were  nearly  all  well  below  our 
average  phase  speed  of  0.51  ms1.  Exceptions  included 
several  isolated  points  for  ease  IV,  which  may  be  the 
result  of  noisy  data  and  near  the  surface  for  case  V;  the 
latter  may  indicate  periods  of  convective  instability. 

Two  Scan  Fish  tows  were  made  in  the  constant  depth 
mode  at  about  30  m  just  east  of  mooring  A2  (Fig.  1h), 
concurrent  with  the  events  of  ease  IV.  A  NLIW  packet 
was  first  crossed  during  a  shoreward-directed  tow  that 
was  immediately  followed  by  a  seaward-directed  tow. 
Shown  in  Fig.  9  is  <rr  from  tows  25  and  26  as  a  function  of 
along-traek  distance.  The  features  marked  a  and  b  cor 
respond  to  case  IV  events  in  the  ADCP  records.  The 
location  of  mooring  A2  is  indicated  by  the  vertical  dot¬ 
ted  line.  The  peaks  at  a  and  b  are  well  defined  by  density, 
and  their  shapes  were  very  consistent  between  the  tows. 
Note  that  the  a  and  b  negative  peaks  arc  about  500  m 
apart  and  are  slightly  closer  together  for  the  seaward 
tow,  w  hich  is  consistent  with  wave  movement  to  the  left. 
Internal  wave  activity  is  apparent  after  the  NLIW  packets 
pass  but  does  not  appear  correlated  between  the  tows. 
Density  data  from  the  string  mooring  (SI;  not  shown) 
show  that  density  contrast  would  have  been  greater  if 
the  level  tow  had  been  made  between  10-  and  20-m 
depths.  Distances  and  corresponding  times  for  the  se¬ 
lected  features  in  Fig.  9  are  listed  in  Table  2.  The  advee- 
tive  correction  here  is  -0.15  m  s  '.  The  data  in  Table  2 
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FIG.  6.  Vertical  velocity  Wat  ADCP  locations  (case  II),  at  lhe  deplh  of  ihe  maximum  VV  time 
variability,  for  time  intervals  showing  three  soliton  events  (labeled  a-c)  propagating  from  A4 
(offshore)  to  A1  (onshore). 

yield  propagation  velocities  of  -0.41  and  -0.43  m  s  1  for 
a  and  b,  respectively.  Correcting  for  the  LF  advcction, 

C  becomes  —0.26  and  -0.28  ms-1.  These  values  are 
consistent  with  ADCP-based  speeds  for  case  IV  but  are 
lower  than  the  averages  in  Fig.  8. 

After  all  moorings  were  recovered,  two  ScanFish  tows 
were  made  near  the  A2  location.  Figure  10  shows  a, 
from  tows  42  and  43  as  a  function  of  along-track  dis¬ 
tance;  we  denote  this  as  case  VI,  which  is  possibly  the 
result  of  a  solitary  wave  breaking  into  a  train  of  dis¬ 
persive  waves.  The  peak-connecting  lines  are  a  best 
guess  at  corresponding  features  in  the  wave  packet,  from 
tow  42  to  43.  As  in  Fig.  9,  the  features  arc  more  com¬ 
pressed  for  the  seaward  tow,  consistent  with  shoreward 
wave  movement.  Furthermore,  the  wave  packet  has 


significantly  evolved  between  tows  42  and  43;  on  the 
other  hand,  it  may  have  interacted  with  another  wave. 
Distances  and  times  (from  similar  plots  versus  time)  for 
the  14  features  yield  14  estimates  of  C.  Note  that  the 
time  intervals  between  feature  positions  in  tows  42  and 
43  increase  shoreward  (with  decreasing  distance)  be¬ 
cause  of  the  respective  tow  directions.  The  C  estimates 
show  an  increase  with  increasing  distance,  but  this  may 
be  biased  by  the  aforementioned  changing  interval. 
The  mean  of  the  14  estimates,  corrected  for  advection, 
is  C  -  -0.54  ±  0.07  m  s"1,  close  to  the  ADCP-derived 
average  speed  of  0.51  ms  '.  MacKinnon  and  Gregg 
(2003)  find  wave  speeds  about  the  same  as  we  find  (about 
-0.5  m  s  l)  from  observations  made  in  19%  near  our 
experimental  area.  However,  Shroyer  et  al.  (2008),  based 
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FIG.  7.  Overall  depth  and  horizontally  averaged  low-frequency, 
rotated  current. 


on  measurements  made  in  2006  in  the  same  general  area, 
report  phase  speeds  for  NLIWs  closer  to  0.7  m  s~l. 

/>.  Wave  properties 

We  now  examine  case  III  in  some  detail,  where  a  pair 
of  solitary  waves  is  seen  moving  past  the  four  ADCP 
moorings.  Horizontal  velocity  U  (rotated  so  that  posi¬ 
tive  U  is  seaward  along  the  mooring  line)  and  vertical 
velocity  W  for  ADCPs  A1-A4  are  shown  in  Fig.  1 1 ;  note 
that  the  maximum  depths  shown  in  these  plots  are  less 
than  the  actual  water  depths  of  67,  65,  76.  and  01  m  for 
A1-A4,  respectively.  The  time  intervals  of  the  plot 
panels  are  0.07  day  (about  100  min)  for  the  respective 
periods  containing  two  propagating  solitons.  The  ob¬ 
served  propagation  speed  of  event  a  is  slightly  greater 
than  that  of  event  b.  Note  that  V  changes  sign  near  20-m 
depth,  near  the  maximum  in  the  magnitude  of  W.  The 
negative  (positive)  U  above  (below)  the  zero  level  is 
consistent  with  propagation  in  the  negative  x  direction 
or  shoreward  along  the  mooring  line  for  a  first-mode 
wave.  Amplitudes  of  the  waves  decrease  at  A2  and  A1 
relative  to  A3  and  A4,  and  smaller  wavelike  features 
appear  between  a  and  b,  suggesting  possible  breaking 
or  fissioning  of  the  solitons  (Djordevic  and  Rcdekopp 
1978).  Temperature,  salinity,  and  density  for  case  III  at 
A2  and  A3  are  shown  in  Fig.  12  (as  stated  above,  max¬ 
imum  depths  shown  arc  less  than  actual  water  depth), 
where  vertical  displacements  of  10-20  m  are  evident, 
particularly  for  temperature  and  density.  The  event 
times  from  the  ADCP  records  are  indicated  in  Fig.  12  by 
vertical  dashed  lines.  The  time  lags  at  the  string  moor¬ 
ings  reflect  the  slightly  different  positions  (—200  m)  of 
the  moorings.  The  appearance  of  two  additional  small 
waves  at  SI  (A2)  is  consistent  with  what  is  seen  in  Fig.  1 1. 

At  the  depths  of  maximum  variance  in  W.  short 
(20  min)  segments  of  data  centered  on  the  points  where 
W crosses  0  (negative  to  positive;  e.g.,  at  /()in  Fig.  6)  were 
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Fig.  S.  Mean  phase  speeds  based  on  arrival  times  of  selected 
internal  soliton  events  in  the  records  of  vertical  velocity  at  the 
ADCP  locations:  (top)  corrected  for  adveetion  and  (bottom)  un 
corrected. 


Meor  phose  speed  (corrected  for  adveetion) 


A4-A3  A3-A2  A2-A1  Overall 

-0.51  ±0.09  m  s 


Mean  phose  speed  (uocorrect e d ) 

A4-A3  A3-A2  A2-A1  Overall 

-0.54  ±0. 1  1  m  s” 


fitted  by  least  squares  to  the  functional  form  for  W  in 
Eq.  (1)  to  determine  the  parameters  a  and  r  (-L/C). 
Recall  [from  Kq.  (1)J  that,  for  Eulerian  data. 
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FIG.  9.  Scan  Fish  level  lows  (top)  25  and  (bottom)  26.  The  location 
of  ADCP  A2  is  indicated  by  a  vertical  dotted  line. 
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Tahlh  2.  'limes  and  distances  for  ease  IV  events  in  ScanFish 
tows  25  and  26. 


low'  25  (a) 

l  ow  26  (a) 

low  25  (b) 

l  ow  26  (b) 

Time  (days) 

176.9945 

177.0294 

176.9966 

177.0312 

Distance  (km) 

13.93 

12.70 

14.44 

13.15 

If  is  the  peak  amplitude  of  the  fitted  W.  then 


awr  / 

a  ~  —-—-/max 


a.  t 


seelrl  - - — \  tanh^ - - 


2  X  0.3849 


(3) 


However,  the  time  scale  resulting  from  the  fitting  algo¬ 
rithm  is  biased  by  the  advection  (as  in  the  determina¬ 
tions  of  C)  and  lienee  needs  to  be  eorreeted  before  (3) 
ean  be  used.  That  is.  if  the  LF  current  is  in  the  direction 
of  propagation,  the  observed  time  scale  r  w  ill  be  shorter 
than  the  intrinsic  r  and  viee  versa.  We  ean  then  write 

Ct  ~  L~  (C  +  U)t  or 

T  =  (l  +  ^V.  (4) 

In  Fig.  13.  we  show  the  20-min  segments  of  W  for  the  two 
events  of  ease  111  at  the  four  ADCP  locations.  On  each 
subpanel  is  listed  the  correlation  of  the  fit  of  (2)  to 
the  data  and  r\  The  amplitudes  elearly  decrease  going 
from  A4  to  Al.  The  least  squares  fitting  procedure  was 
performed  on  data  for  the  five  eases  yielding  a  total  of 
36  determinations  of  aiX  and  r  .  Correlations  for  all  36 
events  are  >0.81  and  most  are  >0.90.  Adveetion  cor¬ 
rections  (4)  were  made  by  using  the  average  value  of 
C  ~  0.51  m  s  1  and  the  displacement  amplitude  a  eal- 
eulalcd  using  (3).  Calculated  r  and  a  are  relatively  in¬ 
sensitive  to  the  variability  of  C  about  the  average  used. 
On  the  other  hand,  L  depends  directly  on  C,  so  the 
variability  of  L  depends  on  that  of  both  C  and  r.  Over¬ 
all  average  soliton  period  (T  -  3.6t)  is  ~10  min,  and 
wavelength  (for  this  C;  A  -  3.6/.)  is  —300  m  (as  in 
Holloway  1987).  Means  and  standard  deviations  of 
a  and  at  each  location  are  shown  in  Fig.  14,  and  the 
means  are  listed  in  Table  3;  as  indicated  in  Table  1,  there 
are  10  samples  each  for  A1-A3  but  only  6  for  A4.  Al¬ 
though  the  variability  at  eaeh  location  is  high,  a  decrease 
in  amplitude  of  both  vertical  displacement  and  velocity 
with  deereasing  water  depth  is  suggested.  This  ehange 
with  depth  eannot  be  directly  explained  with  the  simple 
model  presented  here;  below,  we  examine  bottom  fric¬ 
tion  and  the  potential  for  shear  instability  in  the  pye- 
noeline  to  explain  the  amplitude  decrease.  Furthermore, 


D'stance  (km) 


FlG.  10.  ScanFish  level  lows  (top)  42  and  (bottom)  43.  The  lo¬ 
cation  of  ADCP  A2  is  indicated  by  a  vertical  dotted  line.  Corre¬ 
sponding  NL1W  peaks  arc  connected  by  solid  lines. 

these  displacement  amplitude  magnitudes  \a\  are  al  or 
above  a  maximum  amskX  (means  also  listed  in  Table  3) 
for  KdV  (and  extended  or  second-order  KdV)  waves 
(Helfrieh  and  Melville  2(X)6:  Small  2001). 

Horizontal  velocity  U  suggests  a  two-layer  structure 
(Fig.  1 1),  where  £  is  the  height  of  the  interface  between 
the  layers,  of  thickness  hx  and  h2 .  Here,  £  can  be  esti¬ 
mated  by  integrating  vertical  velocity:  that  is,  £  = 
j\V(-h |,  t)  dt.  For  ease  III,  £  excursions  at  A2  and  A3 
are  comparable  to  the  pycnocline  depressions  in  the 
bottom  panels  of  Fig.  12.  To  simplify  comparison  to 
KdV  theory,  wc  compute  layer-averaged  velocities  U \ 
(upper  layer)  and  U2  (lower  layer).  Beeause  the  zero- 
crossing  depth  varies  among  the  events,  the  depth 
interval  between  15  and  30  m  is  excluded  in  the  calcu¬ 
lation.  Correlations  between  U2  and  £  for  all  36  selected 
events  are  strongly  negative,  consistent  with  shoreward 
propagation.  The  resulting  series  are  fit  to  seeh  [(/  —  /())/r] 
to  determine  amplitude  for  eaeh  event  period.  Average 
amplitudes  for  U\  and  U2  al  eaeh  location  are  plotted 
in  Fig.  I4e.  The  U\  layer  averages  are  quite  noisy  in  the 
latter  part  of  the  record,  and  correlation  values  for  the 
least  squares  fit  are  correspondingly  low.  Hence,  only 
correlations  above  0.8  are  used  in  the  averages  show  n 
in  Fig.  14e.  This  value  represents  the  lowest  correlation 
found  for  the  W  fits,  and  the  U2  correlations  all  exceed 
this  value.  This  problem  is  most  noticeable  for  A4,  where 
only  four  events  have  acceptable  correlations.  The  num¬ 
bers  of  events  included  in  the  averages  are  indicated  in 
the  figure.  The  signs  of  the  averages  show  n  are  consistent 
with  shoreward  movement. 

c.  Energetics 

It  is  useful  to  calculate  energies  associated  with  the  ob¬ 
served  NLIW  events.  Wc  assume  I  hat  the  wave  available 
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Fig.  1 1 .  Rotated,  (left)  high-pass-filtered  horizontal  velocity  £/and  (right)  vertical  velocity  W  for  case  III  from  (top  to 
bottom)  the  four  ADCP  records.  The  features  marked  a  and  b  can  be  tracked  through  all  four  moorings. 


potential  energy  (APE)  and  kinetic  energy  (KE)  are 
approximately  equal;  hence,  the  total  energy  per  unit 
distance  along  the  crest  (perpendicular  to  propagation 
direction)  of  internal  waves  can  he  taken  as  —2  X  APE. 


It  is  likely,  however,  that  the  kinetic  energy  is  greater 
than  the  available  potential  energy  for  large-amplitude 
solilons.  The  ratio  KE/APE  may  be  closer  to  1 .5  (Mourn 
el  al.  2007;  Klymak  and  Mourn  2003;  Klyniak  el  al.  2006), 
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l  ie;.  12.  (top  to  bottom)  Temperature,  salinity,  and  a,  from  conductivity  temperature  string  moorings  (left)  SI  and 
(right)  S2  for  ease  111.  The  time  intervals  displayed  are  as  in  F:ig.  1 1.  Vertical  dashed  lines  mark  corresponding  event 
times  in  ADCP  records  (sec  Fig.  1 1 ). 


so  the  equipartition  assumption  may  underestimate  the 
total  energy.  That  said,  the  total  energy  can  be  expressed 
(in  a  two-layer  system)  as  E.  2APE  -  Apg£2(.v)  dx 

(Phillips  1969),  where  Ap  =  p2  -  p\  is  the  layer  density 
difference  (Table  3)  and  „v  represents  distance  relative 
to  the  moving  waveform.  The  units  of  E{  are  J  m  *.  The 
limits  X\  and  t2  arc  chosen  to  include  the  features  of 
interest;  X\  is  selected  to  be  7  min  (equivalent  distance 
of  210  m  for  a  phase  speed  of  0.5  m  s”1)  prior  to  the  0 
crossing  for  W  (see  Fig.  5),  and  x2  is  selected  at  7  min 
prior  to  the  0  crossing  for  W  in  the  following  event,  to 
include  possible  smaller  waves  resulting  from  breakup 
or  fission  of  the  first  wave  (see  section  3b).  This 


approach  limits  the  calculation  of  energy  to  cases  having 
at  least  two  events  (e.g..  1 1  la— 1 1  lb  yields  one  result  for 
each  of  the  four  ADCP  locations).  In  all,  there  are  five 
energies  each  for  A1-A3  and  three  at  A4.  There  is 
considerable  scatter  in  the  results,  but  averages  for  each 
location  arc  shown  in  Table  3.  Note  that  the  mean  en¬ 
ergies  for  A1  and  A2  are  lower  than  those  for  A3  and 
A4,  suggesting  a  transition  between  A2  and  A3.  Overall 
means  and  standard  deviations  for  these  location  pairs  are 
E\2  =  0.13  ±  0.09  MJ  m  1  and  EM  =  0.31  ±  0.11  MJ  ni 
the  difference  is  significant.  The  NLIW  energies  esti¬ 
mated  here  are  larger  but  of  the  same  order  as  most  of 
those  reported  by  Shroycr  et  al.  (2010)  for  observations 
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FlG.  13.  Hie  unfilled  W  (black)  and  fit  led  (green;  see  text)  at  depth  of  maximum  variance  for  the  two  events  [(top)  a  and  (bottom)  b)  for 
ease  ill  and  (from  left  to  right)  A1  to  A4.  Fitted  amplitude  correlation  with  observed  data,  and  observed  time  scale  r'  are  listed  in 
each  panel. 


made  in  2006  just  to  the  south  of  our  experimental  area. 
The  change  in  energy  might  he  due  to  frictional  effects. 
To  estimate  dissipation  between  these  two  areas,  wc  first 
calculate  energy  flux  averaged  over  the  feature  of  in¬ 
terest,  F  -  CEJ(xz  -  *i)  (Mourn  et  al.  2007).  We  divide 
by  x2  -  -V i  to  give  a  horizontal  mean  versus  integrated 
energy;  Fthen  has  units  of  W  m_I.  The  fluxes  thus  cal¬ 
culated  for  each  event  are  then  averaged  as  above  to 
yield  F]2  =  63.4  W  m  1  and  F?4  =  274.2  W  m  The 
dissipation  rate  is  estimated  as  D  =  (F34  -  Fr)/pl.\h  ~ 
2.75  X  10  W  kg  \  where  Aa  =  10  km  (distance  be¬ 
tween  a  point  midway  between  Al  and  A2  and  a  point 
midway  between  A3  and  A4),  p  —  1025  kg  m  \  and  h  - 
75  m  (the  average  depth  for  the  moorings;  see  Fig.  lc 
for  reference).  This  value  of  D  is  two  orders  of  magni¬ 
tude  larger  than  that  found  by  Hal  lock  and  Field  (2005) 
from  data  in  a  nearby  area,  but  their  value  was  based  on 
a  long-term  mean  flux  and  averaged  over  a  larger  hori¬ 
zontal  area;  they  stated  that  dissipation  rates  averaged 
over  individual  NL1W  events  would  likely  be  much 
greater.  MacKinnon  and  Gregg  (2003)  found  different 
dissipation  values  (about  3  X  10  *  W  kg  ');  they  also 
used  different  spatial  limits  for  averaging.  Shroycr  ct  al. 
(2010)  report  dissipative  loss  values  of 50-100  W  m  1  for 
NLlWs,  comparable  to  our  FM  -  F]2  =  —200  W  m  1 


(they  calculated  dissipative  loss  as  density  times  the 
dissipation  rate  integrated  vertically  over  the  depth  of 
the  pycnocline  and  horizontally  over  the  wave).  The 
dissipation  may  be  the  result  of  a  combination  of  both 
shear  instability  in  the  pycnocline  and  bottom  frictional 
drag  on  the  horizontal  velocity  associated  with  the  waves 
(Shroycr  ct  al.  2010;  Mourn  et  al.  21X37;  Klymak  and 
Mourn  2003;  Dewey  and  Crawford  1088). 

Richardson  number  Ri  calculations  based  on  the  den¬ 
sity  profiles  from  the  string  moorings  (SI ,  S2)  and  velocity 
shear  from  ADCPs  A2  and  A3  show  some  tendency  for 
instability  (Ri  <  0.25)  mainly  shallower  than  25  m  but 
more  at  A2  than  at  A3.  A  quite  crude  estimate  of  dissi¬ 
pation  can  be  made  from  Ri-bascd  mixing  formulations 
(c.g..  Large  ct  al.  1904;  Zaron  and  Mourn  2009).  The 
^-profile  parameterization  (KPP)  formulation  of  Large 
et  al.  (1094)  expresses  the  eddy  diffusivity.  Kp  =  K{)\  1  — 
(Ri/Rio)2]1  for  Ri  <  0.7  and  K(t  =  0  for  Ri  >  0.7,  where 
Rio  ~  0.7  and  Kq  =  5  X  10  3  m:  s  *.  We  noted  that 
diffusivity  is  high  at  SI  compared  to  S2,  suggesting  that 
more  mixing  activities  are  likely  at  SI  than  that  at  S2.  An 
estimate  of  dissipation  rate  can  be  obtained,  because 
Kn  -  Q.lelN2  (Osborn  1080).  Dissipation  rate  is  then  e  ** 
5KnN2<  which  is  on  the  order  of  10  6  W  kg  ',  but  dis¬ 
sipation  is  smaller  for  deeper  depths,  where  Ri  >  I  (for 
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Flo  14.  Mean  (±  I  standard  deviation)  amplitudes  at  each  ADCP  location  (diamonds  con¬ 
nected  by  dotted  lines),  (a)  Displacement  amplitude  a.  (b)  Fitted  vertical  velocity  amplitude  aw. 
The  bathymetric  profile  along  the  mooring  line  is  also  shown  (with  right  ordinate  axis), 
(e)  Upper  and  lower-layer  horizontal  velocity  amplitudes  <?ui  and  y  numbers  in  parentheses 
for  V\  indicate  how  many  events  arc  included  in  the  respective  statistics  (sec  text).  Triangles 
denote  derived  amplitudes  bV\  and  hV2  from  the  two-layer  model  (Table  3). 


most  of  the  water  column),  and  depth-averaged  dissi¬ 
pation  is  more  likely  one  order  of  magnitude  smaller 
than  at  the  shallower  depth  bins.  A  water  column  average 
dissipation  is  then  (assuming  about  5-m-thick  pycnoeline) 
D  shear  -  5e/75  -7  X  10  *  W  kg  \ 

Bottom  dissipation  rates  were  computed  by  assuming 
that  the  “wall  layer1'  or  the  constant  stress  layer  exists 
between  the  sea  bed  and  the  nearest  ADCP  bin  (see, 
e.g.,  Johnson  et  al.  1994;  Wijesekcra  et  al.  2010).  This 
is  a  very  rough  estimate,  because  this  boundary  layer 


thickness  might  be  smaller  than  the  assumed  height 
(i.e.,  to  the  height  of  the  first  bin  of  ADCP,  about  4  m). 
The  friction  velocity  was  estimated  for  a  hydrodynami- 
cally  smooth  bottom  (e.g.,  for  mud  flats)  with  a  bottom 
roughness  length  scale  of  about  Zo  =  0.01  mm  (e.g., 
Sutton  1953).  These  dissipation  estimates  can  be  treated 
as  lower  bounds,  because  z»  can  be  higher  than  0.01  mm. 
Calculated  total  dissipations  exhibit  strong  time  variabil¬ 
ity,  ranging  from  0.1  X  10  6Wkg  '101.2X10  6 W kg  1 
at  —  1 2-h  (and  shorter)  periods,  roughly  corresponding 
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Table  3.  Two-layer  paramclers.  I  he  lop  10  rows  arc  based  on 
observations.  Rows  tl  17  are  derived  using  the  firsl-order  KdV 
model  [Eqs*.  (5)-(8)]. 


A1 

A2 

A3 

A4 

H  (m) 
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lh  (m) 
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«vi  (m  s  ’) 

-0.15 
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0.19 
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0.16 
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0.29 
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Co  (ms  *) 

0.44 
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0.46 

0.44 

C  (m  s  l) 

-0.49 

-0.48 

-0.51 
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A  (m) 
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265 
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T  (min) 

8.97 
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N:i  (ms  ’) 
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/>c:  (ms  ') 

0.10 

0.12 

0.14 

0.12 

E  (MJ  m  ’) 

0.14 

0.18 

0.30 

0.33 

to  semidiurnal  and  internal  wave  periods.  A  similar  cal¬ 
culation  based  on  LF  scries  (tidal  and  longer  periods) 
yields  a  dissipation  of  about  30%  of  the  total.  As  for  the 
thermocline  dissipation,  we  find  an  average  value  on  the 
order  of  10  6  W  kg-1.  Again,  this  value  is  for  the  ~4-m 
assumed  bottom  boundary  layer,  so  it  would  be  about 
10%  of  this  value  if  averaged  over  the  entire  water  col¬ 
umn.  This  yields  about  the  same  value  for  Dk» tt«m  as 
found  above  for  Ds hcar.  The  sum  of  /^bottom  and  D shear  is 
then  about  10  W  kg  \  which  is  lower  than  but  on  the 
same  order  as  the  wave  dissipation  D  computed  above. 

4.  Discussion 

Additional  two-layer  KdV  parameters  for  a  system 
with  layer  thicknesses  h\  and  h2  and  densities  p\  and  p2 , 
respectively,  separated  by  an  interface  £  described  by 
Fq.  (1)  are  (Small  2001;  Holloway  1987;  Gear  and 
Grimshaw  1983), 


k'/i./i,  ,  A  P 

(5) 

W  8  8lh  ' 

(  ah,  ~  h,\ 

C{i+2  V<2  )’  anJ 

(0 

4  h]h\ 

(7) 

3(/j,  -  h2)a ' 

where  a  is  negative  for  solitons  of  depression  with  h\  < 
h2.  Upper-  and  lower-laver  horizontal  velocities  are 

U^C-f  «*“>  "l  =-y£-  <*) 

n\  n2 


with  amplitudes  buy  =  C{)ulh\  and  bV2  =  ~C{)alh2% 
respectively. 

Based  on  the  layer  parameters  (top  four  rows  of  Ta 
ble  3)  and  the  fitted  values  of  (row  5),  the  quantities 
expressed  in  Eqs.  (5)-(8)  were  calculated  and  are  listed 
in  Table  3  (rows  1 1-17).  Note  that  that  magnitude  of  C  is 
about  10%  greater  than  that  of  C<>  for  the  four  locations, 
reflecting  the  nonlinear  contribution  (6)  to  C.  Also,  the 
values  of  C in  Table  3  are  essentially  the  same  as  the  mean 
value  (0.51  ms  l)  found  in  section  3a,  for  the  observa¬ 
tions,  and  do  not  change  significantly  for  the  four  loca¬ 
tions.  The  amplitudes  6Ui  and  b u2  are  plotted  in  Fig.  14e 
(triangles)  for  comparison  with  observational  values: 
these  show  quite  good  agreement  and,  except  for  U2  at 
A4,  fall  within  the  observed  variability.  There  is  gener¬ 
ally  a  slight  increase  in  horizontal  velocity  amplitude 
with  increasing  water  depth,  as  was  found  for  the  ob 
served  values,  but  this  may  not  be  significant.  These  cal¬ 
culations  [(5)— (8)J  are  dependent  on  stratification,  and 
the  CTD  station  associated  with  A4  was  displaced  about 
2  km  northeast  and  conducted  5  days  later.  There  was 
some  evidence  of  a  front  in  the  region,  which  may  add 
uncertainty  to  the  results. 

The  energy  per  unit  crest  length  of  a  two-layer  soliton 
is  given  by  (Holloway  1987;  Osborne  and  Burch  1980) 

l-:  =  ^\f)ga2L.  (l>) 

Calculated  values  of  E  appear  in  Table  3  for  the  four 
locations.  The  E  values  are  comparable  to  those  calcu¬ 
lated  directly  from  the  data  Ec  and  show  a  similar  de¬ 
crease  as  the  solitons  move  into  shallower  water. 

5.  Conclusions 

Six  eases  of  strong  NLI  Ws  were  identified  on  the  outer 
continental  shelf  off  New  Jersey  over  a  cruise  period 
of  about  10  days.  The  phase  speeds  were  affected  by 
the  background  currents,  which  were  primarily  tidally 
driven.  Water  depth  and  stratification  also  affected 
the  overall  characteristics  of  the  waves.  Adveetion- 
eorrected  phase  speeds  averaged  0.51  ±  0.09  m  s=  1  along 
the  mooring  line.  Uneorreeted  phase  speeds  differed 
from  corrected  phase  speeds  by  10%-20%.  Average 
wavelength  was  about  300  m,  and  average  wave  period 
was  about  10  min.  Wave  parameters  were  consistent  be¬ 
tween  the  various  observational  techniques  that  utilized 
both  currents  and  hydrography.  ScanFish  tows  at  a  fixed 
depth  (see  section  2)  provided  better  wave  characteriza¬ 
tions  than  the  profiling  tows,  because  horizontal  spatial 
aliasing  could  not  be  overcome  by  adjusting  profiling 
parameters  and  ship  speed.  Wave  amplitude  decreased 
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in  bolh  vcrlical  displaeemenl  and  velocity  wilh  shal¬ 
lowing  water  dcplh.  Wave  energy  correspondingly  de¬ 
creased  as  lhe  waves  propagated  up  the  sloping  bottom 
into  shallower  water.  The  observed  internal  solitons  can 
be  locally  described  by  an  inviscid,  first -order,  two-layer 
KdV  model,  but  the  apparent  decrease  in  amplitude  as 
lhe  waves  move  into  shallower  water  perhaps  requires 
lhe  use  of  a  second-order  (cubic)  extended  KdV  or 
more  a  complex  model  to  fully  explain  this  decrease.  On 
the  other  hand,  frictional  effects  may  be  important  in 
shallow'  water.  Rough  estimates  of  dissipation  due  to 
bottom  drag  and  shear  instabilities  in  the  pycnocline  are 
of  lhe  same  order  as  lhe  average  dissipation  inferred 
from  the  energy-flux  divergence  near  the  center  of  lhe 
mooring  array. 
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